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Abstract 


The  entropy  of  activation  &^^for  electron  transfer  reactions  occurring 

/V 

In  concentrated  electrolyte  solutions  is  calculated,  taking  into  account  the 
electrostatic  and  electronic  contributions.  By  applying  the  mean  spherical  and 
exponential  approximations  to  treat  a simple  model  of  electrolytes  (hard 
charged  spheres  in  a dielectric  continuum),  the  pair  distribution  function 
is  calculated.  These  approximations  are  also  used  to  calculate  the  elec- 
trolyte reorganization  term,  since  both  quantities  contribute  to  the 

ieltJu  £,  . 

electrostatic  part  of  . Numerical  application  is  made  to  the  ferrous- 

A 

ferric  exchange  reaction  in  0.55  M HC-t04.  Results  of  a preliminary  electron 
tunneling  calculation  for  the  electron  transfer  are  included.  Reasonable 

tfv  is  obtained. 


agreement  between  observed  and  calculated 
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I.  Introduct ion 


It  is  a pleasure  to  dedicate  this  article  to  our  esteemed  and  beloved 
colleague,  Veniamin  Levich,  on  the  occasion  of  his  sixtieth  birthday.  His 
invaluable  contributions  to  electrode  kinetics,  and  his  elucidation  of 
quantum  effects  in  electron  transfer  react  ions,  are  too  well-known  to  all 
of  us  to  need  documenting  here.  It  is  our  hope  that  the  next  few  years  will 

see  him  once  again  able  to  be  active  and  thriving  in  this  field. 

The  aims  in  the  present  paper  are  two-fold:  (1*)  the  investigation  of 

the  large  negative  entropy  of  activation  found  for  electron  exchange  reac- 
tions, such  as  in  Eq.  (1.5)  below,  even  at  appreciable  electrolyte  concen- 
trations, and  (2)  the  introduction  of  recent  statistical  mechanical  methods 
for  electrolytes  into  electron  transfer  theory,  to  treat  salt  effects. 

The  symbols  and  some  definitions  are  first  recalled:  The  rate  con- 

• i 

stant  of  electron  transfer  reactions  is  k ; the  activation  energy  E is 

i 3 

defined  from  the  slope  of  a plot  of  -tn  kr  versus  the  reciprocal  of  the  ab- 
solute temperature  T as  in 

Ea  = - k 9 In  kr/d(l/T)  , (1.1) 

where  k is  the  Boltzmann  constant  (usually  in  cal/mol e/deg).  Eq. 

(1.1)  permits  kr  to  be  written  as 

kr  = A exp( -Ea/kT)  . (1.2) 

A can  thereby  be  calculated  from  the  experimentally  measured  k and  E . 

r a 

The  kr  can  also  be  written  in  absolute  rate  theory  form  as1 

kp  = ( kT/h) exp(  -AG VkT)  = ( kT/h) exp(  -AH*/kT) exp(  ASf/k)  . (1.3) 

A comparison  of  (1.1)  to  (1.3)  then  yields,  using  the  Gibbs -Helmholtz  equation, 
exp(AS^/k)  =A/(kTe/h)  . (1.4) 


AS^  is  the  entropy  of  the  transition  state  minus  that  of  reactants,  each  in 

some  standard  state.  The  value  of  AS^  depends  on  the  units  of  A (e.g.,  A can 

be  In  -t/mole/sec  or  cm3/mole/sec,  say)  and  thereby  the  standard  state  is 
implied  in  AS^  (e.g.,  1 mole/cm3  or  1 mole/i). 
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For  electron  exchange  reactions  such  as 

Fe(H20)6+2  + Fe(H20)6+3  Fe(H20)6+3  + Fe(H20)6+2  (1.5) 

the  experimental  values  of  aS^  are  qui^e  negative,2  e.g.,  -25  cai/mole/ 
deg,  and  E is  9.9  kcai/mole,  when  measured  at  a fairly  large  ionic 

a ~ 

strength  (0.55  M perchloric  acid),  the  excess  perchloric  acid  being  needed 
to  inhibit  hydrolysis  of  the  Fe(H20)6+3  ions.  The  standard  state  chosen  in 
AS^  was  1 M,  by  choosing  the  units  of  A to  be  //mole/sec.  This  large 

negative  AS^  is  quite  common  for  electron  exchange  reactions  between  small 
hydrated  ions3  and  so  our  considerations  are  intended  to  apply  to  other  such 
react i ons  bes ides  (1.5). 

There  are  a number  of  contributions  to  AS^,  including  translational, 
electronic  (any  electronic  nonadiabaticity  produces  a contribution),  and 
electrostatic  in  the  presence  of  the  added  electrolyte.  The  first  of  these 
is  readily  estimated,  the  second  has  been  roughly  estimated  in  a recent 
study,4  and  the  third  is  the  main  concern  of  the  present  paper. 

The  subject  matter  of  the  paper  is  subdivided  as  follows:  In  Section 

II  the  "nonequilibrium"  distribution  of  coordinates  in  the  transition  state 
is  recalled.5  The  various  contributions  to  the  free  energy  of  activation 
AG^  are  given  in  Section  III  and  the  electrostatic  one  is  further  discussed 
in  Section  IV.  In  Section  V,  the  'hiean  spherical"6  and  "exponential"7 
approximations  for  treating  the  electrolyte  effects  are  described,  aqd 
the  equations  used  to  calculate  the  electrolyte  effects  are  summarized. 

An  outline  of  their  derivation  is  given  in  Appendix  A.  The  numerical  results 
are  presented  in  Section  VI  and  discussed  there,  together  with  a comparison 
with  a more  elementary  model. 
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II.  Distribution  Function  in  Transition  State 
The  N-particle  distribution  function  p for  ai  I coordinates  q (transia- 
tions,  vibrations,  and  orientations)  of  all  molecules  has  been  shown9  to 
be  of  the  form  p^,  where 

p*  = C exp(-uVkT)  (2.1) 

C is  normalization  constant  and 

U+(q)  = Ur  + m(Ur-UP)  (2.2) 

in  the  transition  state.  Here,  Ur(q)  is  the  potential  energy  function  of  the  enti 

/v 

system  when  the  reacting  pair  consists  of  reactants,  while  UP  is  that  when  that  pal 
consists  of  products;  m is  a Lagrangian  multiplier,  obtained  from  a free 
energy  balance  equation.  The  argument  leading  to  (2.1)  Is  recalled:  The 

transition  state  lies  at  the  intersection  of  the  U r and  UP  potential  energy 
functions.  To  find  it  one  minimizes  the  free  energy  of  formation  of  this 
state  subject  to  this  constraint.8  The  Helmholtz  free  energy  of  formation 
AF^  consists  of  an  energetic  term  and  an  entropic  term,  minus  that  of  the 
original  reacting  system 

= Jp  ^Urdiq  + kT  J*  p ^ -tn  p ^dq^  - constant  . (2.3) 

The  equations  of  constraint  are 

f P Vdq  = fp  fUPdq  (2.4) 

and 

f Pfdq  = 1 . (2.5) 

Minimizing  (2.3)  subject  to  (2.4)  and  (2.5)  yields 


0 * J*  (Ur+kTtn  p^)  Sp^dq 


(2.6) 
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0 = J (Ur-UP)5pfdq  (2.7) 

0 = J*  fep^dq  . (2.8) 

Multiplying  Eqs.  (2.7)  and  (2.8)  by  Lagrangian  multipliers  m and  -kttn  C 
and  adding  them  to  (2.6)  yields  (2.1)  and  (2.2). 

The  molecular  parameters  appearing  in  Ur  and  Up  include  bond  properties 
(force  constants,  equilibrium  bond  lengths,  anharmonicity  constants),  ionic 
charges,  dipole  moments,  and  Lennard-Jones  parameters,  UrandUP  differ  mainly 
In  the  ionic  charges  of  the  reacting  pair  and  in  the  bond  properties  of 
that  pair.  In  determining  the  distribution  function  p^,  the  charges  of 
the  reacting  pair  behave  as  though  they  were  e,^  and  e2^,  where 

e,^  = e.r  + m(e.r-e.p)  ( i = 1,2),  (2.9) 

ejr  being  the  charge  on  particle  1 when  it  is  a reactant,  and  ejP  being  that 
when  it  is  a product. 

The  above  derivation  used  a canonical  ensemble,  which  Is  related  to 
the  Helmholtz  free  energy  F.  Eqs.  (2.2)  and  (2.9)  also  apply  when  an  iso- 
baric  ensemble,  which  is  related  to  the  Gibbs  free  energy  G,  is  used  in- 
stead. 5C 
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III.  Free  Energy  Changes 

We  consider  first  the  free  energy  of  formation  of  the  transition  state, 
in  which  the  reactants  are  at  some  separation  distance  R.  It  contains  a 
translational  contribution,  AgJ„  d,  an  electronic-nonadiabatlcity  con- 
tribution9  denoted  by  -kT  -In  h,  a work  wr,  both  electrostatic  and  nonelec- 
trostatic,  to  bring  the  reactants  to  a distance  R,  a configurational  con- 
tribution associated  with  forming  the  transition  state  at  a distance  R from 
the  equilibrium  system  of  reactants  at  the  same  R,  AG_^,  and  a term  close 
to  unity  associated  with  the  fact  that  transition  state  may  have  fluctua- 
tions8 in  R,  denoted  by  -kT  -In  P: 

= AGL„5,R  + “r  * V - kT  <•"  Ip  . (3.  I) 

AG^ans  r ar'ses  as  follows:  The  six  translational  degrees  of  free- 

dom of  the  pair  of  reactants  have  a partition  function  ( 2TTm,kT/h2)  3^2 
( 2TTm2kT/h2) 3/2,  where  m,  and  m2  are  masses  of  the  reactants.  These  degrees 
of  freedom  become  three  translations  of  the  transition  state,  with  a parti- 
tion function  ( 2rr[m|-wn2]  kT/h2)  3^2,  two  rotations  associated  with  a moment 
of  inertia  p.R2  [p.  = m,m2/( m,4m2)]  and  a partition  function  8rr2pR2kT/h2, 
and  one  reaction  coordinate  R.  Thereby,  one  finds 

“trans,R  = -kTt"2h/kT  <3-2> 

where  Z equals 

Z = ( frT  kTAi)  i/2R2  . (3.3) 

and  is  the  average  col  I ision  frequencyof  two  particles,  at  unit  concentration, 
moving  in  a continuos  medium.  The  reaction  rate  is  now10 
kr  = pxZexp[-(wr  + AGRfy/kT]  . 


(3.4) 


wr  Itself  can  be  expressed  in  terms  of  the  pair  distribution  function 
g ,2( R) , defined  more  fully  in  Sec.  IV;  g,2  is  unity  (outside  the  hard 
sphere  diameters)  when  the  two  particles  do  not  otherwise  Interact  with  each 
other  or  with  the  medium. 


wr  = - kT  -tn  g ,2(R)  (3.5) 

whence 

k.  = VkZg  , 2( R) exp(  -AGRW)  . ( 3. 6) 

It  has  been  shown  elsewhere  that  when  AG  ^ depends  quadrat ica I ly  on  the 

K 

"fluctuation"  m in  Eq.  (2.2),  it  can  be  written  as5b 

AGRf  = m2(\0+\.)  , (3.7) 

where  m is  the  Lagrangian  multiplier  in  (2.2),  \.  is  the  inner  shell  con- 
tribution (e.g.,  due  to  changes  in  bond  lengths  and  angles  in  the  reactants) 
and  X0  is  the  contribution  arising  from  outside  the  f i rst  coordinat  ion  shel  1 s. 


\o  includes  the  solvent  and  electrolyte  reorganization  outside  those  shells. 

m is  obtained  from  a simple  free  energy  balance  equation.  sa»  5*5 

AgJ-AgJp  =A60l+wP-wr  (3.8) 

R K 

where  ag|P  is  the  same  as  AG^  but  starting  from  the  separated  products,  and 
AG01  is  the  "standard"  free  energy  of  reaction  in  the  prevailing  medium  and 
at  the  prevailing  temperature. 5 

It  has  been  shown  that  the  polar  "outer"  contribution,  m2Xo.  can  be  expressed 
in  terms  of  the  differences  of  free  energy  of  two  equilibrium  systems. 5^,5tl 
Namely,  denoting  m2X0  by  AGR  o*.  The  latter  can  be  written  as 

AGr!o  = m2Xo  = G'-°°P  - G,„0  • (3*9) 

Here,  G|-o  is  the  electrostatic  contribution  to  the  free  energy  of  a system 
whose  central  charges  are  e^-e^l  = 1,2)  and  which  is  otherwise  similar  to 
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the  actual  system,  while  G,_0°P  is  the  corresponding  quantity  for  a system 

|* 

with  central  ion  charges  e.'-e.  but  whose  orientational  coordinates  of  the 
solvent  and  position  coordinates  of  the  Ion  atmosphere  do  not  respond  to 
any  change  in  charges  on  the  two  central  ions;  only  the  electronic  polariza- 
tion of  the  system  responds. 

Since  the  charges  e.^  on  the  central  ions  are  given  by  (2.9),  the 
central  ions  in  system  [ 1 -0]  in  (3.9)  have  charges  -mAe  on  ion  1 and 
mAe  in  ion  2,  where  Ae  is  e,P-e|r,  and  so  is  -(e2P-e2r).  I.e. , 

State  [1-0]  = (-mAe,  mAe)  , (3.10) 

the  charges  being  a distance  R apart. 

To  calculate  the  aG„  in  (3.9)  it  is  necessary  to  calculate  the  elec- 
trostatic  terms  G|-0  °P  and  G , _0 » for  the  above  two  equilibrium  systems. 

For  this  purpose,  as  well  as  for  calculating  wf  in  (3.5),  the  various  con- 
tributions to  the  electrostatic  free  energy  are  listed  in  the  next  section. 

For  the  symmetric  system  (1.5),  treated  in  the  quadratic  approximation, 
m equals  - 1/2,  a result  which  can  be  derived  from  (3.7)-(3.8),  using 
AG° 1 = 0 and  wf  = w^. 


IV.  Electrostatic  Contribution  to  the  Free  Energies  In  Eg.  (3.9) 


As  just  noted,  one  needs  the  free  energy  of  equilibrium  systems. 
Consider  a system  where  the  charges  on  the  two  central  ions  are 
e,  and  e2,  a distance  R apart,  the  charges  of  the  1:1  electrolyte  are 
63  = “®4  = e,  and  the  medium  has  a dielectric  constant  D.  There  are  two 
contributions  to  the  electrostatic  free  energy  G. 


G = GS°,V  + Gatm  . 


(4.1) 


where  GS°^V  is  the  solvent  contribution,  outside  of  the  inner  shell,  defined 
as  the  contribution  in  the  absence  of  the  added  electrolyte.  (The 
inner  shell  contribution  to  GSO*v  is  contained  in  m2\.  in  (3.7)  and  not  in 
the  m2\0  of  (3.9)  and  so  is  not  included  in  (4.1).) 

In  a simple  electrostatic  model,  Gs°'v  is  the  Born  charging  term  plus 
the  worke,e2/0R  required  to  bring  the  ions  together.11 


where  a,  and  a2  are  the  radii  of  the  solvated  ions. 

One  can  see  immediately  how  Eq.  (4.2)  can  be  used  in  (3.9).  One 
replaces  e,  and  e2  by  the  charges  in  system  [1-0],  namely  the  -mAe  and  mAe 
of  Eq.  (3.9),  then  replaces  D by  DQp.  One  subtracts  from  this  result  the 
corresponding  quantity  with  D replaced  by  Dg.  Thereby,  the  solvent  con- 
tribution to  AG^  is 


fsol  v 
R,o 


. !)(_!_ . -L) 

RmD  « ' 


op 


(4.3) 


a well-known  result.  5,12 

We  turn  next  to  the  Gatm  in  Eq.  (4.1).  Gatm can  be  wri tten  as  the  electro- 
static free  energy  of  i nteraction  of  ions  I and  2 with  their  ion  atmospheres  when 


m 


mmm 


1 and  2 are  far  apart,  plus  the  change  in  this  ion-atmosphere  electrostatic 
free  energy  when  the  ions  1 and  2,  surrounded  by  their  atmospheres,  are 
brought  from  an  infinite  separation  distance  to  a distance  R apart.  We 
denote  these  two  contributions  by  Gatm(°°)  and  AGatm(R). 


Gatm  = GatmH  + AGatm(R)  . 


.atm 


(4.4) 


Both  terms  in  G should  reduce  to  zero  at  zero  ionic  strength  by  defini- 
tion, and  so  the  work  term  e,e2/DR  of  (4.2)  is  not  included  in  Gatm. 

To  calculate  Gatm  it  is  convenient  to  define13’ 14  a pair  distribution 
function  g . j( r)  for  particles  i and  j:  If  Pj  is  the  average  number  of 

particles  of  type  j per  unit  volume  then  g . j ( r)  is  defined  by  the  statement 
that  the  average  number  of  particles  of  species  j in  a volume  element  4rrr2dr  at  a 
distance  r and  r+dr  from  a particle  of  species  i is  p^g . j(  r)4TTr2dr.  g.j(r) 
has  the  property  that15 


g.j(r)  -•  1 as  r » . 


(4.o) 


Consider  a charging  process  of  the  central  ions,  where  the  charge  of  ions  1 
and  2 is  Xe,  and  \e2  and  X goes  from  0 to  1.  The  charges  in  the  ion  atmo- 
sphere are  e3  and  e4.  The  g , j when  1 and  2 are  infinitely  far  apart  is 
denoted  by  g,j(r;\e,)  at  stage  \ of  the  charging  process  (j  = 3,4).  The 
potential  acting  on  ion  1 at  its  center,  due  to  the  ion  atmosphere  is, 

= J{  (P3e39  1 3 (r;Xe,)  + p4e4gl4  ( r ;Xe  ,)  ]/Dr}4TTr2dr  . (4.6) 


Multiplying  by  the  charge  increment  d(Xe|),  integrating  over  X,  and  adding 
the  corresponding  term  for  ion  2,  one  has  Gatm(») : 

Gatm(=o)  = / (♦  ,Xe , + f2Xe2)dX  . 

X=0 


(4.7) 
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Finally,  the  term  AGatm(R)  in  (4.4)  can  be  expressed  in  terms  of  ap- 
propriate differences  of  electrostatic  work  terms  to  bring  particles  I and 
2 together  in  the  medium:  The  electrostatic  contribution  of  the  atmosphere 

can  be  defined  as  the  total  work  to  bring  particles  1 and  2 together  when 
they  have  charges  e,  and  e2  minus  this  work  when  the  atmosphere  is  uncharged, 
i.e.  , minus  -kT  -tn  g 1 2 °(  R)  + e|e2/DR.  Using  (3.5)  one  has 

AGatm(  R)  = - kT  -tn[g  1 2(  R)/g  1 2°(  R)]  - (e,e2/DR)  . (4.8) 

At  zero  ionic  strength  this  AGatm(R)  vanishes. 

Since  the  ion  atmosphere  in  the  [1-0, op]  system  in  (3.9)  does  not  res- 
pond to  the  charges  on  the  central  ion,  the  atmosphere  contribution  to  the 
Gi-o0*5  in  (3.9)  is  zero,  and  so  the  atmosphere  contribution  to  AG^  ^ is 
-G|-o,  and  is  the  sum  of  (4.7)  and  (4.8),  with  e,  and  e2  replaced  by  -mAe 
and  mAe,  respectively. 

It  remains  to  calculate  the  9jj's  'n  (3.6),  (4.6)  and  (4.8),  and  for 
this  purpose  the  mean  spherical  and  exponential  approximations  are  described 
in  the  next  section. 

Before  giving  these  equations,  it  is  useful  to  just  consider  a more 
elementary  model,  based  on  Debye-Huckel  theory.  The  solution  of  the  linear- 
ized Poisson-Bol tzmann  equation  for  the  total  electrostatic  potential  at  a 
point  r outside  of  a central  ion  1,  fi(r),  is16 

'h(r)  = ( e |/Dr)exp[-x(  r-R)]/(  1+xR)  ( r > R)  ( DH)  , (4.9) 

when  R is  the  distance  of  closest  approach  of  the  atmosphere  to  the  central 
ion.  DH  denotes  Debye-Huckel.  The  electrostatic  work  required  to  bring  an 
ion  2 of  charge  e2  to  a distance  R,  W)2(R),  is  e2i|i,(R),  i.e., 


_ 


II 


W,2(R)  = (e,e2/DR)/(  1+hR)  ( DH)  . (4.10) 

The  value  for  the  pair  distribution  function  g)2(R)  is,  in  the  Debye-Huckel 
approximation, 

g.t(R)  = 1 - PWI2(R)  , (4.11) 

where  p = 1/kT.  If  one  used  the  exact  W|2(R),  its  relation  with  gi2(R) 

would  be 

9 l2(  R)  = exp(-pWl2(R))  (4.12) 

If  one  were  to  use  (4.10)  and  (4.11)  one  would  encounter  a serious  error-- 
even  a negative  g 1 2 ( R) . If,  °n  the  other  hand,  one  were  to  use  (4.10)  and 
(4.12),  i.e. , an  "exponential"  form  of  the  Debye-Huckel,  the  error  in  g,2 
is  found  to  be  less.  17  In  treatments  of  kinetic  salt  effects  by  Debye- 
Huckel  theory,  one  in  fact  forgets  about  the  linearization  (4.11)  used  there 
in  the  derivation,  by  using  (4.12)  and  focusing  attention  on  the  calcula- 
tion of  a free  energy  term  and  including  it  in  the  free  energy  of  acti- 
vation, i.e,,  including  it  in  an  exponential  in  the  rate  constant  expression. 

An  uncertainty  arises  in  the  Debye-Huckel  picture  when  the  diameter 
of  the  principal  ions  in  the  atmosphere  differs  from  that  of  central  ion 
2.  In  the  interest  of  simplicity  this  difference  in  diameters  will  be 
neglected  in  Eqs.  ( 4-1 3) -( 4. 1 7) , just  as  it  was  in  (4.10). 

In  computing  the  -tn  g 1 2/g  12°  term  in  (4.8),  where  now  g 1 2 0 is  unity, 

one  obta i ns 


-kT  In  gl2(R)/g  |2°(R)  =WI2(R)  =(e  ,e2/DR)/(  1+xR)  ( DH) , (4.13) 

regardless  of  whether  one  uses  (4.11)  and  then  expands  the  logarithm,  thus 
"undoing"  that  approximation,  or  uses  (4.12)  directly. 

The  Ga*m(°°)  term  given  by  (4.7)  and  (4.6)  is  also  estimated:  The  gl3 

appearing  in  (4.7)  is  given  by 
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g,3  = J - Pe34f ,(  r)  = 1 - e30\e  ,{exp[-x.(  r-R)]  ]/(  l4*R)Dr  . (4.14) 

Noting  that  g u is  given  by  a similar  expression  with  e3  replaced  by  e4, 
i.e.,  by  -e,  Eqs.  (4.6)-(4.7)  yield 

Gatm(a>)  = -(e,z4«22)K/2D(  1 +kR)  . (4.15) 

Finally,  one  obtains  from  (4.4),  (4.8),  (4.13)  and  (4.15),  that 

Gatm  = - x[^e,24e22)+e,e2]/D(  1+kR)  ( DH)  . (4.16) 

The  reorganization  term  is  obtained  by  replacing  e,  by  -mAe  and  e2  by 
mAe,  whence 

AGR,oatm  = 0 0>H)  • 

Other  Debye-Huckei  methods5*3  give  slightly  different  results  than  (4.15)- 
(4.16),  e.g.,  from  solution  of  the  linearized  Poi sson-Bol tzmann  equation 
when  ions  1 and  2 are  fixed  a distance  R apart.  Gatm  is  then  defined  as 
the  electrostatic  work  to  charge  these  two  ions  minus  that  when  their  ion 
atmosphere  is  absent. 


V.  Mean  Spherical  and  Exponential  Approximations 


A.  Description  of  Approximations 

The  mean  spherical  approximation  ( MSA)6  employs  the  linearization  pre- 
sent in  Debye-Huckel  theory  and  introduces  hard  sphere  radii,  but  is  sta- 
tistical mechanical  in  nature.  Whereas  the  Debye-Huckel  theory 
does  not  allow  for  hard  sphere  repulsion  between  the  ions  belonging  to  the 
ion  atmosphere,  the  MSA  does.  The  MSA  yields  good  agreement18  with  activity  co- 
efficients, calculated  numerically  (Monte  Carlo),  in  the  concentration 
range  0.1  tol  M.  It  does  not  give  accurate  results  for  the  radial  distri- 
bution functions  on  contact.70' 1 8 For  this  purpose,  the  "exponential  approxima- 
tion" ( EA) 1 which  makes  use  of  the  MSA  result  but  attempts  to  correct  for 
the  linearization,  is  much  better,  as  judged  by  comparison  with  exact  statis- 
tical mechanical  numerical  (Monte  Carlo)  computations.70 

To  calculate,  thereby,  the  g,2(R)  at  contact  in  (3.6)  and  (4.8)  the  EA 
is  used  below,  while  to  calculate  (4,6)-(4.7)  the  MSA  Is  used.  (For 
activity  coefficients  its  results  are  comparable  with  EA.) 

These  approximations  are  described  next. 

First,  recalling  the  definition  of  g.j(r)  in  the  previous  section  one 
notes  that  when  there  is  no  interaction  between  the  particles  then  for  all 
i and  j 

9 1 j ( 0 =1  (all  r)  (5.  1) 

Thereby,  h j j ( r)  defined  as 

h.j(r)  = 9.j(r)  - 1 (5.2) 

reflects  the  interparticle  interactions  and,  in  virtue  of  (4.5),  approaches 
zero  for  large  interparticle  separations. 
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It  Is  very  convenient  to  define  a new  funct  ion,'  s Cjj(  r) , related  to 
hjj(r),  but  that  expresses  correlations  between  i and  j more  directly 
than  does  g.j,  by  being  more  related  to  the  intermoiecuiar  potential  it- 
self. For  that  one  decomposes  h(  r) , for  a one  component  fluid,  into  a 
chain  of  integrals  of  this  new  function  C(  r) . That  is,  if  C(rl2)  correlates 
directly  the  small  volume  elements  from  positions  1,  2,  where  r^  denotes  f r • — r 2 f , 
then  h(rt2)  is  constructed  from  the  sum  of  the  direct  correlation  C(rl2) 
plus  the  direct  correlation  of  a particle  at  r,  with  the  particles  p dr3 
in  r3,  multiplied  by  the  direct  correlation  of  those  at  r3  with  those  in  r2, 
integrated  over  all  r3,  plus  the  direct  correlation  of  the  particles  at  r, 
with  those  at  r3,  multiplied  by  that  of  those  at  r3  with  those  at  r4,  and 
then  by  those  of  r4  with  those  of  r2,  integrated  over  all  r3  and  r4 , and 
so  on. 

h(  r 1 2 ) = C(  r |2)  + Jc(  r 1 3)  PC(  r32)dr3  + J*C(  r 1 3)  pC(  rj4)  pC(  r42)dr3dr4  + ... 

(5.3) 

or  using  the  symbol  * for  a convolution  integral: 

h(r)  = C(r)  + C*PC  + OpC*PC  + ...  (5.4) 

This  definition  of  C( r)  can  be  seen  to  be  completely  equivalent  to 

h(  r)  = C(r)  + h(r)*pC(r)  . (5.5) 

which  gives  (5.4)  upon  successive  iterations.  Eq.  (5.5)  is  the  cele- 
brated Ornstein-Zerni ke  equation.19 

The  generalization  for  mixtures  is  given  as 

hij<rij>  =cij<riJ)  + J pJcii(rit)hij(rtj»‘'Trtj!<lrij  • <5-6> 


In  matrix  form  one  has 
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where  P is  a diagonal  matrix. 

When  each  molecule  is  treated  as  having  a hard  spherical  core,  the 

g..  vanishes  when  r..  is  less  than  the  mean  diameter,  (d.+d  )/2. 
3ij  ij  ' J 

Hence 


1 

2 


(d.+dj.) 


(5.  8) 


Ifone  had  a system  consisting  of  a mixture  of  uncharged  hard  spheres,  the 

interaction  between  molecules  is  zero  outside  the  distances  of  closest 
approach.  Making  the  approximation  that 

C.j(r)  =0  for  r -»  d.^  (hard  spheres)  (5.9) 

and  using  condition  (5.8)  one  obtains  the  well-known  Percus-Yevick  (PY)  ap- 
proximation, 20  which  is  an  excellent  description  for  the  fluid  region  of  hard 

sphere  systems.  This  fluid  serves  as  a reference  system  for  realistic  dense 
fluids  and  for  electrolytes. 21 

To  solve  the  set  of  equations  (5.7)  for  electrolytes  some 
knowledge  of  C.j(r.j)  is  needed  for  those  systems,  just  as  it  was  in  the 
PY  case,  a knowledge  hidden  in  the  higher-order  distribution  functions.  We 
proceed  to  motivate  and  give  the  choice  made  in  the  case  of  the  mean  spheri- 
cal approximation. 

In  the  limit22  of  £ -»  0,  Eq.  (5.4)  shows  that  C approaches  Ih.  But  h 
approaches  the  simple  Boltzmann  value  when  g -•  0: 


-tim  C. .(  r . .)  =h..(r..)  =exp(-(3V,.) 
p_0  'J  'J7  'Jv  1 J7  W 


(5.10) 


where  p denotes  1/kT  and  V.j  is  the  interaction  potential  of  i with  j. 
For  sufficiently  low  (3V.j,  expansion  of  the  exponential  yields23 


Cjj(  rjj)  - • 3V-j(  r.j)  . (5.11) 

The  mean  spherical  approximation,  introduced  by  Lebowitz  and  Percus, 
consists  of  assuming  Eq.  (5.11)  for  r > d . j , for  any  densi  ty  p,  then  solving  (5.6) 
subject  to  (5.8).  The  solution  reduces  to  the  linearized  Debye-Huckel 


Jfe." 
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result  when  d.j  Is  made  to  approach  zero,  and  to  the  hard  sphere  PY  equa- 
tion when  the  charges  e.  are  made  to  vanish. 63 

The  exponential  approximation  of  Chandler  and  Andersen7 attempts  to 

undo  the  linearization  embodied  in  (5. li)  Its  distribution  function 
EA  . 

9 1 j 's  obtained  to  be  the  hard  sphere  pair  distribution  function  9jj° 
multiplied  by  a suitable  exponential,24 


91JEA(rij)  =9Ij°(riJ)exp(-Sw1J) 


(5.12) 


By  an  optimization  procedure  they  obtained  an  effective  potential  W.^, 
which  in  our  case  is  well  approximated  by  the  mean  spherical  approximation 
to  g.j  and  by  g.j0: 

-pWjj  = 9ijMSA  " 9ij°  • (5.13) 


(In  Eqs.  (4.10)  and  (4.12)  of  Sec.  IV,  Eq.  ( 5. 1 2) -(5. 13)  were  also  used 

MCA 

for  1-1,  j=2,  at  rl2=R,  but  with  g,  2 replaced  by  the  Debye-Huckel 

value  (4.11)  for  gi  2,  and  with  g? 2 replaced  by  unity.) 

Some  qualitative  i ns i ght  i nto  (5. 12) -( 5. 13)  isobtained  as  follows: 

9 ? j( r)  *s  related,  via  an  exponential,  to  the  work  required  to  bring  i 

and  j to  a separation  distance  r from®.  Eq.  (5.12)  presumes  that  this 

work  is  the  sum  of  the  hard  sphere  work  plus  an  additional  ionic  work. 
MSA 

Remembering  that  g.j  is  a linearized  g.j,  °ne  would  write  for  it,  in 
the  same  spirit  as  (5.12), 


9ljMSA'  WBWlj  - 9|J°  - 0W(J 


(5.14) 


where  we  have  set  g}j. 0 equal  to  unity  in  the  second  small  term  on  the  ri^it  hand  side. 
However,  (5.14)  is  seen  to  be  none  other  than  (5.13).  Eqs.  (5.12)  and 
(5.13)  provide  the  "exponential  approximation"  and  yield  an  improved 

value  for  the  value  of  g^.  at  r.j  = d j ^ . That  is,  to  use  this  "exponential 
approx i mat iorf1  it  suffices  to  solve  the  problem  in  the  "mean  spherical  approxi- 
mation." 


B.  Equations  for  the  g^'s 

The  g ,2°(R)  for  the  hard  spheres  ( PY)  system,  In  the  presence  of  an 
uncharged  atmosphere  of  particles  3 and  4,  occurs  in (4.8)  and  is  given  by29 

9.2°(R)  = (I-?)-'  +(3R/2L)§(1-?)-2  , (5.15) 

where  5 is  the  ratio  of  the  total  hard  sphere  volume  of  the  atmospheric 
ions  4tt(L/2)3N/3  to  the  total  volume  V: 


? = ttPL3/6  , p = p,  + p4  . (5.16) 

The  W 1 2 present  in  Eq.  (5.12)  and  (5.13)  at  r,2  = R,  Is  shown  by 
Eq.  (A. 28)  and  (A.  30)  of  Appendix  A to  be  given  by 

W,2(R)  = (e,e2/DR)(l+p)2/q2  , (5.17) 

where 

p = ( 1+2  hL) 1 ^2  , q = 1 + p + kR,  h = (4npe2/0kT) 1 ^2  , (5. 18) 

K being  the  inverse  Debye  length.  Clearly  (5.17)  tends  to  e|e2/DR  when 
the  atmospheric  number  density  p tends  to  zero.  The  g , 2( R)  in  Eq.  (3.6) 
is  now  given  by  (5.12),  using  (5.15)  and  (5.17). 

In  Eqs.  (4.6)  a quantity  such  as  9 is(  r) -g  u(  r)  appears  when  es  = 

-e4  = e.  According  to  two  of  the  equations  In  (A. 6),  and  the  relation 
between  g.^  and  h.^  in  (5.2),  this  g 1 s'-g  i«  is  given  by  the  following,  where 
the  charges  on  the  central  ions  are  explicitly  indicated  as  Xe,  and  Xe2, 

9 is(  r;Xe ,)  -gu(r;\e,)  =2h,(r;\e,)  (5.19) 

and  similarly 

92s(  r;Xe2)  - g24(r;Xe2)  =2\e2h,(r;Xe,)/Xe,  . (5.20) 
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(The  latter  could  also  have  been  written  as  2h((r;\e2),  since  h|  is  a 
linear  function  of  the  charge  of  the  central  ions.)  Further,  the  integral 
of  hi  over  r,  multiplied  by  various  constants,  Is  the  same  as  the  integrals 
of  the  gij-gi*  and  g23-g24  that  appear  In  Eqs.  (4.6)-(4.7),  and  is  given  by 
Eq.  (A. 19).  In  this  way  one  obtains  from  (4.7) 

GatmH  = -(e,2+e2z)H/Dq  , (5.21) 

a result  which  may  be  compared  with  (4.15), 

On  the  other  hand  the  AGatm(R)  in  (4.8)  is  obtained  from  the  contact 
gf2  values.  By  using  (4.8),  (A. 28),  and  (A.  29),  one  finds 

AGatn\R)  = (e,e2/DR)[(l+p)2-q2]/q2  . (5.22) 

Thus,  from  ( 4.4),  (5.21)  and  (5.22)  one  has 

Gatm={-(el24e22)(HR/q)+  e ,e2  [(  1+p)  2-q2]/q2}/DR  . (5.23) 

To  obtain  the  atmospheric  contribution  to  AG_  * one  uses  (3.9)  and 
(5.23),  but  replacing  the  et  and  e2  in  the  latter  by  -mAe  and  mAe,  and 
the  D by  D$.  This  gives  the  amosphere  contribution  to  G ,_<>•  There  is 
no  atmosphere  contribution  to  G°^0,  since  by  definition  the  ion  atmos- 
phere and  orientational  coordinates  do  not  rearrange  in  the  [1-0, op]  system; 
only  the  electronic  polarization  responds  to  changes  in  charges  on  the 
two  central  ions.  Thereby,  the  atmospheric  contribution  equals 

AG*a*m-  -Gat_m0  = m2(Ae)2(HR)2/q20sR  . (5.  24) 

In  summary,  AG*  is  obtained  by  adding  (5.24)  to  (4.3),  and  the 
R,  o 

g,2(R)  In  (3.6)  Is  given  by  (5.12),  (5.15)  and  (5.17),  while  Debye-H(/ckel  values 
are  given  by  (4.13)  and  (4.12).  The  former  equations  are  just  as  easy  to 
apply  as  the  latter. 
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Finally,  the  wr  appearing  in  Eq.  (3.5) 
via  (5. 12). 

wr  = W ,2(  R)  - kT  In  g |j°(  R)  . 


is  related  to  V,t(R)  and  g,a°(R) 

(5.25) 


The  expressions  obtained  from  these  equations  for  various  contribu- 
tions to  AS^,  using  -SAG^/dT,  are  as  follows: 

The  translational  contribution  is  obtained  by  differentiating  Eq. 
(3.2),  using  (3.3),  and  is  given  by 

AStrans,R  = k ^Zh/kT)  - ( k/2)  . (5.26) 

The  electronic  contribution  is 


AS 


el 


f 


= k In  k 


(5.27) 


neglecting  any  temperature  dependence  of  x. 

The  entropic  work  term,  -dwr/dT,  obtained  by  differentiating  (5.25) 
and  (5. 17)  is 


AS  * = k In  g |2°  - ( e ,e2/D_TR)[(l+p)/q2]  {0(  1+p) -[xR(  0-l)(  l+p-*L)/pq]  } , (5. 

w 

where  0 denotes  -d  -tn  D /d  In  T. 

s 

The  corresponding  term  at  zero  ionic  strength  is  obtained  by  setting 
p = q/2  = 1.  The  corresponding  Debye-Huckel  term,  obtained  by  different- 
iating (4.10),  is 

AS  r^(  DH)  = -( e ,e2/D  TR)  ( 20+0xR-WtR)/2(  1+xR)  2 . (5.29). 

w 

The  solvent  reorganization  term,  obtained  by  differentiating  (4.3),  is 

ASR,oSO,V  = m2(Ae)20/DsTR  (5.30) 

neglecting  the  temperature  dependence  of  D . 
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The  ion  atmosphere  reorganization  term,  obtained  by  differentiating 
(5.24),  is 

^sJ,oatm  = ”m2(^e)2((HR)2/q3DsTRHqe+(0-i)(i+p-(xL/p)]]  . (5.31) 

The  ion  atmosphere  reorganization  term  in  the  Debye-Huckel  approxima- 
tion, obtained  from(4.i7),  is  zero. 

For  the  electron  exchange  reaction  (1.5),  the  value  of  m is,  as  al- 
ready noted,  -1/2. 
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VI.  Numerical  Calculations  and  Discussion 

Using  the  ionic  diameters26  given  in  Table  I the  various  contribu- 
tions to  AS^  are  calculated  from  the  AG^  contributions  using  -^AG^/dT. 

They  are  summarized  in  Table  I.  The  agreement  with  the  net  experimental 
value  at  0.55  M HC-t04  is  seen  to  be  quite  reasonable,  based  however  on  the 
rough  estimate  for  the  nonadiabat icity  factor  h in  ref.  4.  There  is  seen 
to  be  a substantial  contribution  of  the  ion  atmosphere  to  AS^,  as  expected 
from  the  extensive  atmosphere  shielding  involved.  The  various  solvent  and 
atmosphere  terms  are  further  classified  into  "work"  type  (wr  type)  and  re- 
organization type  (AG^  type)  contributions.  The  g 1 2 °( R)  appearing  in 
Eq.  (5.15)  is  1.08  at  0.55  M HC-tO*. 

It  is  interesting  to  compare  the  effect  of  the  ion  atmosphere  on  AG^ 
and  AS^,  calculated  in  Table  I,  with  that  calculated  from  the  Debye-Huckel 
theory.  The  results  are  given  in  Table  I and  are  seen  to  be  strikingly 
close  to  the  mean  spherical  results  (in  the  exponential  approximation  for 

both).  Nevertheless,  the  MSA  is  readily  generalized  to  more  complex  sys- 
tems, whereas  the  DH  is  not. 

In  the  present  treatment,  the  central  ions  were  taken  to  be  of  the  same 
size  (diameter  R) , as  were  the  ions  in  the  atmosphere  (diameter  L).  One 
can  readily  generalize  to  arbitrary  sizes  of  each  of  these  ions,  using  the 
exponential  approximation,  with  the  form  of  the  mean  spherical  approximation 
results  obtained  by  Blum  and  Htfye. 27^ 

A second  generalization  Is  in  calculating  the  fluctuations  In  separa- 
tion distance  R in  the  transition  state,  instead  of  largely  restricting  It 
to  the  "hard  sphere  contact  distance."  One  would  maximize  with  respect  to 


r,  x(  r)Z(  r)g  t2(  r)exp( -pAG^  ),  r being  the  separation  distance  (replacing 

K,  O 

R).  9 1 2 ( •")  is  given  approximately2715  by 

9 12(f)  =gi2°(r)exp{  H 1 2(  R)( R/r) exp[ -2r(  r-R) ] ] , (6.1) 

where  H|2(R)  is  given  by  Eq.  (A. 29)  of  Appendix  A and  where  2r  ?s  p-1. 

Eq.  ( 6.  1 ) applies  to  equal  size  central  ions  and  a symmetric  atmosphere, 
but  can  be  generalized.27^5 

The  present  equations  become  less  accurate,  as  judged  by  comparison 
of  the  exponential  approximation  with  Monte  Carlo  results, in  the  range 
around  0.01  M for  very  highly  charged  signs.  18  The  phenomenon  usually 
termed  "ion-pairing"  occurs  there,  and  at  lower  densities,  and  one  needs  a 
more  sophisticated  approximation,  such  as  perhaps  the  hyper netted-cha in 
one.  1 8^ 

It  should  be  emphasized  of  course  that  the  model  employed  here  for 
salt  effects  is  the  primitive  one,  containing  as  it  does  only  hard  spheres, 
without  specific  ionic  interactions  and  not  allowing  for  solvent  structure 
breaking  and  forming  effects.  Some  discussion  of  deviations  from  the  di- 
electric continuum  approximation  is  given  in  Ref.  18g. 

In  summary,  the  use  of  the  Debye-Huckel  treatment  to  calculate  the 
electrostatic  free  energy  tacitly  makes  use  of  the  "exponential  approxima- 
tion" but  with  a Debye-Huckel  gt2  replacing  a mean  spherical  one.  The 
mean  spherical  equations  for  symmetrical  systems  are  not  more  difficult  to 
use  than  the  Debye-HUckel  ones.  They  treat  the  hard  sphere  excluded  volumes 
in  a more  self-consistent  manner,  and  can  be  generalized  to  less  symnetrl- 
cal  systems.  The  calculated  result  for  the  entropy  of  activation  in  Table 
I agrees  well  with  the  experimental  one,  but  the  electronic-nonadiabat ic 
contribution*  should  be  regarded  as  a preliminary  one. 
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Appendix  A.  Mean  Spherical  Approximation.  Physical  and 
Mathematical  Ideas. 

A simple  model  of  electrolytes  of  N species  will  be  used,  character- 
ized by  the  number  density  of  ions  p.,  ion  charge  e.,  hard  core  diameter 
d.  and  dielectric  constant  D.  There  is  also  the  added  condition  of  over- 
all electrical  neutrality,  i.e. , 

N 

E Pkek  = 0 . (A.  I) 

k=1  K 

The  M.S.A.  consists  of  Ornstei n-Zerni ke  equation19  plus  the  Lebowitz- 
Percus  assumpt ion6*3  about  the  behavior  of  C.j(r),  namely  Eqs.  (5.7),  (5.8), 
(5.11),  and  Eq.  (A.  1),  with  V.j  being 

V. . = e.e./Dr..  , r..  -»  d..  . (A.  2) 

ij  ' J U U 'J 

One  can,  incidentally,  show  directly  from  these  equations  that  the  M.S.A. 
satisfies  the  St i 1 1 inger-Lovett2 8> 29  moment  relations,  a result  which  fol- 
lows from  an  observation  of  Groeneveld,  30  and  which  has  the  implication  of  pre- 
dicting oscillations  in  the  charge  density  around  a given  ion  at  high  ionic 
strengths. 

The  solution  of  the  equations  (5.7),  (5.8),  (5.11),  (A.  1),  and  (A.  2) 
was  first  obtained  by  Waisman  and  Lebowitz  and  can  be  found  in  references 
6c,  6d  and  6e.  The  techniques  employed  extended  Werthelm's 
method3  1 for  solving  the  PY  equation  for  pure  hard  spheres,  later  used  by  Lebowitz25 
for  the  PY  mixture  of  uncharged  hard  spheres.  For  the  case  of  symmetric 
electrolytes,  the  problem  was  solved  completely,  while  for  the  asymmetric 
case,  the  solution  led  to  a set  of  coupled  algebraic  equations.  Blum, 27 
utilizing  Baxter's32  application  to  fluids  of  the  Wlener-Hopf  Integral 
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equation  solving  technique,  was  able  to  reduce  the  problem  to  solving  a 
single  algebraic  equation  of  higher  order,  an  equation  which  reduces  to  a 
quadratic  one  in  the  case  of  the  present  problem  of  a symmetric  electro- 
lyte surrounding  the  two  central  equal-size  ions.  For  this  simpler  prob- 
lem we  shall  use  instead  a method  of  solving  the  equations  which  avoids 
the  complications  of  the  previous  methods  and  uses  the  known  solution  for 
symmetric  electrolytes  found  in  Ref.  6d. 

In  the  present  case  the  ionic  atmosphere  consists  of  a symmetric  elec- 
trolyte, composed  of  ions  labelled  (3,4),  with 

e3  = -e4  = e , d3  = d4  = L < R , p3  = p*  = p/2  (A.  3) 

and  with  the  central  ions  labelled  (1,2)  with  charges  e,  and  e2,  and 

d,  = d2  = R , p,=p2=0  . ( A.  4) 

Because  of  this  symmetry  it  readily  follows  from  the  Ornstei n-Zerni ke 
equations,  plus  the  associated  conditions,  that  the  matrices  h and  £ are 
symmetric  and  that  for  the  symmetrical  1:1  electrolyte  one  also  has 

1^3  3 = » C33  = C44  . (A.  5) 

These  quantities  are  the  same  as  those  in  the  absence  of  ions  1 and  2, 
since  pt  = p2  = 0.  They  can  be  taken  as  known,  since  the  solution  to  the 
problem  of  a 1:1  symmetric  electrolyte  by  itself  is  available, 33  and  is  used 
later  via  the  functions  a and  h in  (A. 11),  etc.  One  next  defines 

h = (h33-h34)/2  , h = (h33+h34)/2  , C = (C33-C34)/2  , C = (C33+C34)/2 


It  can  be  seen  that  by  adding  and  subtracting  the  Ornstei n-Zerni ke 
equations  for  the  indices  13,  14,  23,  24,  using  (5,11),  (A. 5)  and  the  sym- 
metric nature  of  the  matrices  h and  C,  one  obtains. 


h 1 3 = h|  + h,  , h23  = h|  + (e2h,/el) 


hu  - fi|  - h|  , h2*  - h|  - (e2hj/e,)  , 


and  an  identical  set  of  equations  with  the  h's  replaced  by  C's.  Here, 


the  functions  ht,  C|  and  h|,  Ct  satisfy  the  following  decoupled  linear 


integral  equations: 


h | — C i + P h ,*C 


h|=0,  r < R | ; C,  = -gee, /Dr  , r > R, 


A A A A 


h | — C | + p h ,*C 


h,  = -1  , r < R,  ; C,  = 0 , r > R,  , 


where 


R . s -A  L+R)  . 


Equation  (A. 8)  contains  no  ionic  charges  and  actually  is  the  same 


as  the  PY  approximation  for  a binary  mixture  of  uncharged  hard  spheres 


of  two  species,  one  of  diameter  L at  a density  p and  one  of  diameter  R 


at  a vanishing  density;  h,  represents  the  unlike  particles  total  pair 


correlation.  Its  solution  is  obtained  as  the  limit  of  the  Lebowitz 


general  solution25  of  the  PY  approximation  for  mixtures,  by  letting  one 
of  the  p's  vanish,  but  is  not  needed  here,  and  so  is  not  cited. 


It  remains  to  solve  Eq.  (A.  7).  Auxiliary  functions  CT|(r)  and  cr(  r) 


are  first  defined: 


C | — gee  ,/Dr  r < R , 

a((r)/2npr  = < 


a(  r)/2rrPr  = 


-C  - 0e2/Dr 


r < L 


r > L 


(A. II) 


(Auxiliary  functions  have  been  introduced  to  solve  other  integral  equations, 
and  were  used  for  the  first  time  by  Wertheim31  to  solve  the  PY  equation. 

The  addition  of  the  Coulombic  term  was  one  of  the  key  steps  introduced  by 
Waisman  and  Lebowitz6C"e  to  solve  the  mean  s pherical  integral  equations 
for  electrolytes.)  Using  bipolar  coordinates  Eq.  (A. 7)  can  then  be  written 


6i(r)  = - 


oo  L 2 oo  r+y 

J*  or ,( y) 0(  L-f  r-y | ) J*  o(x)dx  - J*  0|(y)dy  j dx  , (A.  12) 


where 


=1  if  z>0  , 0=0  otherwise  . 


(A.  13) 


One  now  takes  the  Laplace  transform  of  (A. 12),  obtaining 


where 


G i(  s)  = - [(x2se,/2e)  - s2F ,(  s)  - x2V  |]/[s2-x2-sF2(  s)]  , (A. 14) 


Gi(s)  = J a,(  r)e'srdr  ; F,(s)  = j a i(  r)e's  rdr  , 


(A.  15) 


Ms)  = Ks)  -F(-s)  with  F(s)  = J a(r)e~srdr  , 


(A.  16) 


V)  = G ,( 0)  = J a,(  r)dr  . 


(A.  17)  * 


Eq.  (A.  13),  which  contains  two  unknowns  Ft(s)  and  G|(s),  is  a solvable 
functional  equation.34  By  imposing  the  appropriate  analytic  properties35 
of  G|(s)  and  F|(s),  both  of  these  functions  can  be  obtained.  The  details 
are  found  In  the  1 iterature.  6C> 69  By  now  performing  the  Inverse  Laplace 
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transformation,  one  obtains  Oi(r)  and,  thereby  from(A.IO),  the  following 
explicit  result  for  Ct(r): 


2pe2V  j/D  , 


r<^HX 


C,(r)  ={  (20e2V,/D)[l+{V(r-\)V2r}]  , \ < r < R, 


(A.  18) 


■0e  ,e/Dr  , 


r > R i 


where  the  quantities  V and  V|  are  associated  with  the  mean  potential  energy 


of  the  Coulomb  potential  and  are  given  by:36 
00 

V = J 2TTprh(  r)dr  = -x/(  I + p + xL) 
L 


(A. 19) 


V|  = J 2rrprh|(r)dr  = - e,x/eq 


(A.  20) 


where  x,  p and  q are  defined  in  (5.16)  and  (5.18). 

The  explicit  knowledge  of  C,(r),  given  (A.  1 8) -(A. 20) , can  be  used  to 
obtain  h|(r)  by  Fourier  transforming  the  top  line  of  Eq.  (A. 7).  One  ob- 
tains 

K.(k)  = C,(k)/[l-pC(k)]  , (A.  21) 

where  h,(k)  stands  for  the  integral  of  h ,(  r)exp(  ik-r^ r2sin0d0dcpdr,  inte- 
grated over  0,  cp  and  r,  and  similarly  for  C(  k)  and  Ct(k).  Also  one  can 
write  another  important  quantity,  by  noting  that  directly  from  (A. 7)  it 
follows  that  h|(r)-C|(r)  is  continuous  across  r = R ( , yielding  h|(R|+)  » 

C i ( R i +)  - C |(  R | ),  because  h t(  R t ) = 0.  (R,+  = Cim  R,+c  as  e -*  0 for 


e > 0. ) One  thus  obtains 

h|(Ri+)  = -20ee,/DR  ,q 


(A.  22) 


The  component  h,2  can  be  written  now  in  terms  of  known  quantities: 
h 1 2 = H 1 2 0 + H \ 2 (A. 


2 


and  similarly  hu  and  h)2  can  be  expressed  in  terms  of  H|2°  and  Ht2.  Here, 
H|j°  is  defined  by  (A. 24)  and  H|2  by  (A. 26). 37 

H |2  0 = cp°  + ph  i*C  i 

H 1 2 0 = - I , r < R (A. 24) 

9°  = 0 , r ^ R . 

So  defined,  H|2°  satisfies  the  PY  equation  for  the  uncharged  hard  sphere 
mixtures  problem,  but  for  the  central  ion-central  ion  correlation,  and  its 
solution  Is  completely  known. 37  In  particular, 

9 i2°(R)  =1  + Hi2°(R)  = ( 1-S)"'  +|  (R/L)?(  l-?)-2  , (A. 25) 

a result  quoted  in  (5.15). 

Hi2  is  defined  via 

H 1 2 = V + ( e2/ e |)  ph  |*C  | 

(A.  26) 

H|2=0,  r < R ; cp  = - pe,e2/Dr  for  r > R . 

The  main  problem  is  to  obtain  from  Eq.  (A. 26),  using  (A. 22) -(A. 24) , 
the  value  of  H|2(r)  as  r approaches  R from  above.  This  can  be  done  in  a 
cumbersome  direct  way,  or  indirectly  by  first  showing  that  in  the  interval 
(0  < r < R)  h|*C|  is  a linear  function  of  r.  Since  cp(  r)  equals  -(e2/e,)h|*C 
for  r < R,  one  then  finds  a value  for  cp(r):38 

CP(  r)  = (pee2/D)(2V,  + (req,2/e,x2)  , r<R  , (A.27) 

and 

q,  = 2rrpRlh,(R,+)  . (A.  28) 


^■aassflii  ’■ 


From  this  one  obtains 


H,2(R+)  =cp(R+)  - <¥<R")  =-(pe,e2/DR)(I+p)l/q2  . 

From  the  first  equation  of  (A. 23),  expressed  in  terms  of  the 
stead  of  the  h's  and  using  the  first  half  of  (A. 25),  one  has 

Hij(R+)  = 9 i *(  R)  - 9,2°(R) 

Thus,  the  -0W,2  in  Eq.  (5.14)  is  given  by 

- 3W,2(R)  = H,2(R+) 

at  r = R.  Eqs.  (A.29)  and  (A.31)  yield  (5.17). 


(A.2< 
•s  In- 
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26.  (a)  A.  F.  Wells,  "Structural  Inorganic  Chemistry"  (Oxford  University 
Press,  New  York,  1950),  also  ref.  4;  ( b)  R.  H.  Stokes  and  R.  A. 
Robinson,  J.  Soln.  Chem.  2,  173  (1973);  R.  Trfolo,  J.  R.  Grtgera,  and 
L.  Blum,  J.  Phys.  Chem.  80,  1858  (1976);  (c)  B.  B.  Owen,  R.  C.  Miller, 
C.  E.  Milner  and  H.  L.  Cogan,  J.  Phys.  Chem.  65,  2065  (1961). 

27.  (a)  L.  Blum,  Mol.  Phys.  30,  1529  ( 1975);  (b)  L.  Blum  and  J.  Hdye  (to 
be  publ ished). 

28.  F.  H.  Still  Inger  and  R.  Lovett,  J.  Chem.  Phys.  48.  3858  (1968);  49. 

1991  (1968).  ' “ ’ 

29.  If  Qj(r),  a charge  cloud  density  around  an  Ion  of  species  I at  the 
origin,  is  used  to  denote  4n*  £ p.e.g..(r),  then  the  two  moment  re- 


* • 


1 at  ions  are  that  ( i)  J*Q. .( r) dr  from  r = 0 to  *>  is  -e.,  and  that  (li) 


J\(r)r*dr,  multiplied  by  p^e^,  integrated  over  r and  summed  over  -C, 
equals  -6  £ p^e^2/*2  where  k is  the  Inverse  Debye  length.  Stllllnger 

and  Lovett2 8 have  given  strong  physical  arguments  for  these  relations, 
and  the  relations  agree  indirectly  with  Monte  Carlo  results  [cf. 

Ref.  18  and  J.  Rasaiah,  Chem.  Phys.  Letters  7,  260  (1970)]. 

The  moment  relations  predict  charge  oscillations,  as  one  can  show  from 
standard  inequalities.  E.g.,  for  the  symmetric  case  of  a binary  mix- 
ture, p,  = p2  = p,  e,  = -e2  = e,  d,  = d2  = d,  and  the  moment  relations 
read  jQ,(r)dr  = -e,  JQ. i(  r) r 2dr  = -6e/x2,  the  integrations  being  from 
r = d to  r = od,  and  one  finds  from  those  results  that  Q ,( r)  is  not  al- 
ways of  one  sign,  i.e. , charge  oscillations  occur,  when  xd  > JZ  . 


30. 


The  two  moment  relations  are  satisfied  in  the  linearized  Debye-HUckel 
system  with  d.j  =0,  but  only  the  first  is  satisfied  by  the  Debye-HUckel 

theory  with  finite  d.^. 


Groeneveld  (private  communication  to  Waisman)  pointed  out  that  the 
moment  relations  are  satisfied  by  any  approximation  whose  Cjj(r) 


behaves  as  -pe.e./Dsr  at  sufficiently  large  r.  Thereby,  they  are 
satisfied  by  the  Debye-HUckel  equation  when  d.^  = 0.  (When  d}j  4 0 
the  "modified"  Debye-HUckel  does  not  satisfy  the  Ornstein-Zernike 


equations.)  The  moment  relations,  and  the  resulting  charge  oscilla- 
tions do  not  impose  a severe  restraint  on  the  thermodynamic  properties, 
since  the  latter  are  reasonably  well  satisfied  by  equations  such  as  the 
01 ivares-McQuarr ie  one,**  which  do  not  show  the  oscillations.  The  dif- 
ferences might  affect  other  properties,  however,  perhaps  such  as  the 
capacitance  behavior  of  an  electrode  immersed  in  an  electrolyte. 


31.  cf.  M.  Wertheim,  Phys.  Rev.  Letters  10, 321  (1963);  cf.  E.  Thiele,  J. 

Chem.  Phys.  39,  474  (1963). 

32.  R.  J.  Baxter,  Aust.  J.  Phys.  21_,  563  (1968);  R.  J.  Baxter,  J.  Chem. 
Phys.  52,  4559  (1970);  B.  Nob  1 e , "Met hods  on  the  Wiener-Hopf  Technique," 

(Per gamon  Press,  New  York,  1958). 


' 


33.  The  solution  is  given  by  Eqs.  (2)  and  (3)  of  reference  6d,  by  noting, 
however,  that  we  are  using  indices  3 and  4 instead  of  the  1 and  2 
for  the  symmetrical  electrolyte  in  the  atmosphere. 

34.  Actually  in  our  case  is  quite  easier  than  the  general  case,  because 
F2(s)  is  a known  function  and  (A. 12)  and  therefore  (A. 14)  are  just 
scalar  equations  and  not  matrix  ones. 


35.  G,(s)  Is  regular  for  Re  s ^ 0 and  F,(s)  is  entire  in  the  s-complex 

plane. 
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36.  This  Is  a general  characteristic  of  the  M.S.A.  : the  excess  potential 

energy  enters  as  a coefficient  In  the  algebraic  equations  one  obtains 
as  described  before. 

37.  since  ht  and  Ct  are  kncwn,  the  first  and  third  equations  in  (A. 24) 
give  H 1 2 0 for  r > R,  while  the  first  and  second  equations  give  <p°  for 
r < R.  Thereby,  H|2°  (using  the  second  equation)  is  known  for  ail  r, 
and  so  is  <P°.  Eq.  (A.  25)  is  obtained  from  Lebowitz1  solution, 25  as  the 
limiting  case  of  the  pair  distribution  function  of  the  larger  particle 
with  another  of  its  kind  at  the  limit  of  vanishing  density  for  them, 
but  at  a given  finite  density  for  the  smaller  spheres  species, 

38.  H i 2 is  kn<*m  once  h|  and  C|  are  known,  following  a similar  reasoning 
to  that  In  Ref.  37.  Also,  since  cp(  r)  is  a linear  function  of  r for 

0 < r < R,  it  is  completely  determined  by  cp( 0)  and  ^£(0);  l.e. , in 

virtue  of  Eq.  (A. 26),  by  the  values  of  -h,*C,  and  - ^h,*C,)  at  the 

origin.  These  derivatives  are  easily  found  in  terms  of  the  quanti- 
ties given  in  (A. 20)  and  (A. 22). 
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Table  I 

Ferrous-Ferric  Electron  Exchange  Reactions. 
Contributions  to  AS*  and  Related  Data.* 


Contribution 

AS* 

(0.55  M HC-tO*) 
( cal/mole/deg) 

a£ 

(0  M HC£04) 

( cal/mole/deg) 

Translation 

-10. 1 

same 

Nonadiabatic 

-9.2 

same 

Solvent  reorganization 

0.6 

same 

Entropic  work  term  (-0wr/dT) 

-5.2 

-14.9 

Atmosphere  reorganization 

-0.  1 

0 

Total 

-24.0 

-33.6 

Experimental 

-25.0 

— 

Entropic  work  term,  Debye-Huckel, 
(-dwr/aT) 

-5.0 

-14.9 

Atmosphere  reorganization  ( DH) 

0 

0 

AG* 

AG* 

(0.55  M HC-CO*) 

(0  M H0t04) 

( kcal/mole) 

( kcal/mole) 

Work  term  wr 

1.2 

3.2 

Work  term  wr,  Debye-Huckel 

1.2 

3.2 

Atmosphere  reorganization 

0.  02 
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